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be constant. The large L linear asymptotic behaviour of —AE{L) sets in only after a quite 
large retardation time. The suppression of —AE(L) is partly due to the fact that gluon 
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medium compared to vacuum. The latter radiation spectrum is sensitive to the plasmon 
modes of the quark-gluon plasma and has a rich angular structure. 
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1. Introduction 



Jet quenching has long ago been suggested as a possible signal for the quark-gluon plasma 
(QGP) and this has triggered a considerable activity both on the experimental and 
theoretical sides. The theoretical determination of parton energy loss has been the subject 
of many studies, and numerous observables sensitive to jet quenching are currently analyzed 
in the Relativistic Heavy Ion Collider (RHIC) experiments. 

In a (static) thermal or dense medium of very large size L, a parton of high (but finite) 
energy E undergoes a radiative energy loss AE ra d which increases as a power of the energy 
(AE ra d oc LE in the Bethe-Heitler limit and AE rac [ oc L\[E when the non-abelian LPM 
effect is at work |^]), whereas its collisional energy loss AE co n behaves at most logarithmi- 
cally 1 , AE co n oc LlogE J!], |3|, H|. For a large size medium we thus expect AE ra ^ to be the 
dominant source of energy loss for energetic partons. 

In the opposite limit, namely for a parton of asymptotic energy E — > oo crossing a 
medium of finite size L < L cr (but still thick enough, L/X S> 1, where A is the parton 
mean free path), AE rac [ becomes independent of E, AE rac [ cx E°L 2 H, Rj. The finite size 
limit L < L cr should be relevant in practice since a simple numerical estimate within the 
model of Ref. |§ gives L cr = J XE/m 2 D ~ 7 fm for E = 10 GeV, where m^ 1 is the Debye 



screening length in the medium. It was also shown H that finite opacity (L/X > 1) and 
finite E effects lead to a substantial suppression (together with some energy dependence) 
of AE rac i, when compared to the asymptotic ^-independent result. We can thus expect 
AE co u to compete with AE ra d in the case of a finite size medium. Indeed, recent studies 
g, 0, [T§ suggest that for 'jets' of energy on the order of 10 GeV such as those measured 
in AA collisions at RHIC, the collisional energy loss might be comparable to the radiative 
one, both for light [10] and heavy [11] partons. 

Thus an accurate determination of AE co u for large E and finite L is needed in order 
to interpret the suppression of inclusive large pt hadron production at moderate energies 



observed at RHIC in AA collisions Jl2| , 13 1. In previous studies of parton collisional energy 
loss in a QGP, the loss 2 —AE(L) suffered by the hard parton travelling the distance L 
in the plasma was assumed to be given by (—dE/dx)^ times L, with (—dE/dx)^ the 
rate of energy loss per unit distance occurring in a stationary regime, i.e., long after the 
energetic parton has been produced. This is certainly a valid approximation when the 
medium size becomes very large, and formulas for (—dE/dx)^ of a heavy quark of mass 
M>T travelling in a QGP of high temperature T have been obtained [j|, f|]. 

In the present paper we argue that such an approximation to —AE(L) is incorrect 
in the case of an energetic parton produced at an initial time t = inside a thermal or 
dense medium of moderate size L, which is the relevant situation when discussing jet or 
hadron quenching at large px- As in previous studies of collisional energy loss [^, ||, we 
work in the theoretical high temperature small coupling limit j < 1 for the QGP (and 
in the hard thermal loop - HTL - resummation framework [14, and assume a fixed 



lr The logi? factor arises when AE co ii is calculated in the fixed coupling approximation Q B, 0, but is 
expected to be absent in the case of a running coupling M. This point is briefly discussed in section 3.2. 



From now on all energy losses will be implicitly collisional, — AE = —AE co 
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coupling a s . We also consider a static (non-expanding) QGP in thermal equilibrium. 
As argued in section 3.2, we expect the retardation of the stationary regime found in 
this paper to be qualitatively unchanged in the case of a running coupling. However, a 
rigorous treatment taking into account the running of a s would be needed in order to 
obtain a reliable quantitative estimate of —AE(L). Recalling moreover that g(mrj) > 1 in 
realistic phenomenological applications, we stress that our results should be considered on 
a qualitative level only. 

For a fast parton prepared at t = — oo and travelling in an infinite medium, the 
collisional energy loss can be understood as follows. The proper (chromo-) electric field of 
the parton polarizes the medium, which creates an effective (medium-induced) electric field 
around the parton, responsible for its slowing down. If the energetic parton is produced at 
t = (via some process involving a hard scale ~ pt), we may expect the rate of collisional 
energy loss —dE/dx to be reduced during the time the parton proper field is created, before 
reaching the asymptotic value (—dE/dx)^. A result suggestive of such a retardation effect 
is qualitatively what we find in the following. A fast parton produced initially in the 
medium needs to travel some distance before losing energy at the highest rate. Our main 
conclusion is that collisional losses used in the analysis of nuclear modification factors Raa 
at large pr (for jet or hadron production) have been systematically overestimated due to 
the neglect of this effect. 

In section 2 we present our model for the induced collisional energy loss —AE(L) of 
a parton produced initially in a QGP. In section 3 we give our numerical results for the 
L-dependence of —AE(L), which display a strong attenuation of the energy loss when 
compared to previous estimates, as well as a large retardation time t re t of the stationary 
regime. We give a simple interpretation of the largeness of i re t- We also discuss the domain 
of validity of our analysis, and how the running of a s could affect —AE(L). Finally, in 
section 4, we show that within our (standard) definition of collisional energy loss, —AE(L) 
gets a contribution from induced gluon radiation which arises from the sudden acceleration 
of the parton at t = 0. The in-medium bremsstrahlung spectrum due to charge acceleration 
is indeed not the same as in vacuum since it is sensitive to the longitudinal and transverse 
plasmon modes of the QGP. The angular spectrum presents a diffraction pattern depending 
on the plasma size L. The induced radiated energy is negative, and thus contributes (but 
only partly) to the suppression of — AE(L) and to the retardation effect discussed in section 
3. We conclude in section 5. 



2. Model for collisional energy loss of a parton produced at t = in a QGP 

We derive in this section the master equation ( |2.S| ) for the collisional energy loss —AE(L) of 
an energetic (and massive) parton of velocity v, initially produced in a quark- gluon plasma, 
and travelling the distance L in the medium. We first present our model for the classical 
partonic current density = (p a ,j a ), and then give, in the abelian approximation, the 
expression of the (chromo-) electric field induced by this current density. 
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2.1 Model for partonic current 

In previous studies of parton collisional energy loss in a QGP ||, the classical color 
charge has been assumed to be produced at t = — oo. In this case the current 4-vector in 
coordinate space reads: 

j£?(t,S)=q a V»6 3 (x-vt) , (2.1) 

where denotes the parton 4- velocity, V = (l,v)- The color index a is carried by the 
parton color charge q a defined by q a q a = Cro s , where Cr is the Casimir invariant of the 
color representation R the parton belongs to (Cr = Cir = 4/3fora quark and Cr = Ca = 3 
for a gluon). In 4-momentum space K = (uj, k) the current (|2.1|) becomes 



j^(K) = 2irq a V»5(K.V) , (2.2) 

which trivially satisfies current conservation, K.j^ = 0. 

When considering the quenching of large pt jets or hadrons in relativistic heavy ion 
collisions, the large pt parent parton is rather created at t = in the underlying hard 
partonic subprocess. Its associated classical 'current' is thus, instead of (2.1), 



ti a (t,x)=q a V^(x-vt)0(t) , (2.3) 

giving in momentum space: 



VP 

j^(K)=iq a = q a V» 

K.V + irj 



iP[ YV )+7T5(K.V) 



(2.4) 



However, the 'current' (2.4) is not conserved, K.jq ^ 0. Color charge conservation requires 
the hard parton to be produced in conjunction with at least another parton in the partonic 
subprocess. As a generic conserved partonic current, we will consider the simple case of 
a (color singlet) dipole produced at t = 0, consisting of two partons (with the same color 
charge q a ) of 4-velocities V\ = (1, «i) and = (1,^2): 

f a = { a f) = ig a 1 _ 2 (25) 

\K.Vi + it] K.V2 + ir]J 

In the following we will consider the second parton to be a static heavy quark, i.e. v*2 = 0, in 
which case the space component of the current is simply given by the first term of ( |2.5| ) . The 
choice (^^) instead of ( p.4|) is however crucial. Current conservation indeed constrains the 
form of the electric field (^^) , and allows to treat unambiguously the potential singularity 
at uj = (see ( |2.6| ) and ( |2.7D ), as we briefly explain below. We expect the main results of 
our study not to depend on the simplifying assumption vi = 0. 

2.2 Induced electric field 

Following 0] (see also |16|]), in the abelian approximation and within linear response theory 
the Maxwell equations can be solved in 4-momentum space K = (uj, k), yielding the total 
(chromo-)electric field E a in terms of the classical vector current density j a : 

e L E a L + (e T - k 2 lu?)E a T = —(f L +j a T ) . (2.6) 

iuj 
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The longitudinal and transverse components are given by jl = (j.k/k 2 )k (we denote k = 
\k\) and jt = j — 3l- in the linear response approximation the longitudinal and transverse 
dielectric functions of the plasma and €t are not affected by the external current. We 
consider a high temperature QGP, for which €l.t have been obtained in |17], 18] and later 
rederived in the gauge-invariant HTL resummation framework [jlj, |15[ . 

In the following we will have to deal with the 1/ui potential singularity appearing in 
Q2.60 , and specifically affecting the longitudinal part of the electric field. The latter actually 
arises from Coulomb's law k €lEl oc p, where p is the charge density, by using the equation 
for current conservation p = kji/ui. More precisely, our conserved current (2.5) satisfies 
p = kji/{uj + ir/), for any v\ and V2 (including V2 = 0). This shows that the potential 
1/ui singularity appearing in (^1]) should be regularized with the retarded prescription, 
ui — > ijj + irj. 

We obtain from (I2.6T) the medium-induced electric field £ a , 



£ a {t,x) 



00 du 

-oo w 



d 3 k 
4n 3 i 



-i(ujt—k.x) 



ACL ao, 
■>L_ _| JT 



(2.7) 



Lad 



where j a is given by ( |2.5[ ) (with = 0), or equivalently by ( |2.4D (with v = vi), and the 1/uj 
singularity should be treated with the retarded prescription. In ( |2.7D the subscript denotes 
the implicit subtraction of the vacuum contribution (corresponding to = ey = 1). Since 
the dielectric functions (and the external current) are real quantities in coordinate space, 
implying in momentum space 6l,t(—K) = €l,t{K)* (and a similar relation for the current), 
the expression fl2,7|) is easily checked to be also real. 



2.3 Master equation for parton collisional energy loss 

During its travel in the plasma between t = and t = L/v, the induced energy gain AE of 

the parton of constant velocity v\ = v equals the work of the electric force on its trajectory, 
o-r r L / v 



namely AE = v ■ f Q dt q a £ a (t, x = vt) or: 



AE(L) = q a v- 



00 duj 



d 3 k f L ' v 

47T 3 i 



dt e 



-iK.Vt 



Aa -a 
J_L_ _| JT 

e L e T - k 2 /ui 2 



(2.8) 



ind 



We stress that the latter expression is valid in the abelian approximation for the hard 
parton dynamics and within linear response theory, implying that | AE(L)\ should be small 
compared to the initial parton energy E, which is also consistent with the assumption of a 
constant velocity. 

the current (^.2|) is used in (2.8), the exponential factor equals unity, and 



If as in 

the collisional energy loss on the distance L is uniquely determined by its rate per unit 
distance (—dE/dx)^ = —AE/L. Inserting instead (|2.5|) in ([2.S|), and performing the time 
integral, we obtain the collisional energy loss —AE(L) of a hard parton produced at t = 
and travelling the distance L in the medium, 



We note that the factor between the curly brackets in (|2.9| ) can be rewritten as 



Using the following identities, 
sin(uL) 



(K.vy 



ir5(u) 



K.V 



sin 2 (uL) 



K.V 



nS(u) , 



(2.10) 



(2.11) 



we find that in the L —* oo limit, the expression ( pTTOj) is equivalent to 2n S(K.V)L/v, 
and (|2.9|) thus reproduces the result for (—dE/dx)^ obtained in Ref. ||. In the small L 
limit, however, the second term of fl2.10| ) is subleading and ( [2.10D reduces to (iP(l/K.V) + 
tt6(K.V)) L/v as can be seen also directly from Q2.9Q , thus leading to a modification of 
collisional energy loss at finite L. 



3. The retardation effect 

In this section we first express the energy loss Q2.9D in terms of the discontinuity (on the 
real axis) of the longitudinal and transverse thermal gluon propagators. We then present 
and discuss our numerical results, which show that the asymptotic large L behaviour of 
—AE(L) sets in only after some retardation time. 

3.1 Expression of —AE{L) in terms of thermal gluon spectral densities 

The dielectric functions can be expressed in terms of the longitudinal and transverse ther- 
mal gluon self-energies 3 , 



e L = l + U L (x)/k 2 ; e T = l-Il T {x)/LU 2 , 
where x = to/k and Hl,t have been obtained in the HTL approximation [14, 15 

'x + 1 



(3.1) 



U L (x) = m D 



1 



x i 

o log , 
2 v x 



1 



n T (a:) 



12 2 

2in D x 



1 



2x 



1, fx + 1 

- los {— 



(3.2) 



The Debye mass denoted as mo is given by m 2 D = 47ra s T 2 (l + nj/6) (with nj = 2 the 
number of thermally equilibrated flavours). We also use the longitudinal and transverse 
thermal gluon propagators 

-1 -1 



Al(u> = kx, k) 



At(w = kx, k) 



k 2 + U L (x) ' 1K ' ; LO 2 - k 2 - U T (x) 

Using (3.1) and ( |3.3| ) the expression (2.9) becomes (with v 2 = v 2 cos 2 9, Vj, = v 2 sin 2 9): 
-AE(L) 



(3.3) 



CrO! s 



-IV 



j4 / — [k 2 cos 2 9A L (u;,k)+uj 2 sm 2 9A T {uj,k)}. nd 



■y g— i(ui— kv cos 9) L/i> 



(cj — A;i> cos 0) (lo — kv cos + irf) 



(3.4) 



^In the following, we will use the sign conventions and notations of Ref. . 
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The longitudinal and transverse thermal gluon propagators have singularities on the 
real w-axis, namely branch cuts (corresponding to Landau damping) in the spacelike \x\ < 1 
region, and poles corresponding to collective excitations of the plasma (plasmons) in the 
timelike |x| > 1 region. Those singularities must be treated using the retarded prescription 
uj — > uj + irj arising in the analytical continuation from imaginary to real frequencies in finite 



temperature field theory. As we explained in section 2.2, the potential 1/uj singularity at 
uj — > must be also regularized with the retarded prescription. Thus all singularities on 
the real axis appearing in ( |3.4j) are implicitly written with the same, retarded prescription. 

It is convenient to perform the w-integral in ( |3.4[ ) using Cauchy's theorem, by closing 
the integration contour in the lower (complex uj) half-plane, as required by the presence 
of the exponential factor in the integrand (since L > 0). Thus the integral over the real 
axis is replaced by the integral over the (clockwise-going) contours around the singularities 
which lie just below the real axis, namely the poles at uj = —irj and uj = kvcos9 — irj, and 
the plasmon poles and cuts of the propagators Al,t- We obtain from fl3.4|) : 

-AE{L) . 2 f d 3 k ( 1 



C R a s 



iv 2 j ^ {(-2»7r)l(l - e ikLcose ) Re A L (0, k) 



COS $ 

+(-2z7r) k \ReA L (kvcos9,k) + v 2 sin 2 9 Re A T (kv cos 6, k)] 

/oo 
duj [k 2 cos 2 9lmA L (uj,k) + uj 2 sin 2 8 Im Ay(w, k)] 

/l\ / 1 \ 1 _ p -i(u-kv cos 9) L/v 1 

x p (- p I . (3.5) 

\uj ) \uj — kvcosvj uj — kv cos 6 

) ind 

Due to the fact that Re Al^{uj, k) is an even function of uj, the second line of the latter 
equation vanishes after angular integration. Recalling |[^] that the spectral densities (with 
s = L or T) 

p s (uJ, k) = 2 Im A s (uj + irj, k) = 2ir sgn(oj) z s (k)8(uj 2 - uj 2 (k)) + /3 s {uj, k)0{k 2 - uj 2 ) (3.6) 
vanish at uj = 0, we can rewrite ( |3.5| ) in the form: 

—AE(L) f d 3 k ( I- cos (kL cos 6) 2 [°° duj 



/dk ( 1 — cos (kL cos 9) 9 f°° duj r , 9 9 ^ 9 0n n 

2^{ +v [ fc2cos2 ^^- 2 - 2 ^-] 



x 2 sin 2 ((^-^cosg)L/(27;)) | 

(w - A;v cos 6) 2 J ind ' 



The expression (|3.7] ) for the collisional energy loss is actually ultraviolet divergent. 
The logarithmic UV divergence ~ j dk/k appears in the (leading) asymptotic L — > oo 
behaviour of (|377|), 



-AE^L) _ y 2 f d 3 k f°° duj 



Cro> s 



l d ^IZ^ [k2cos29pL+uj2WepT] ^ 



irL 

— 5(uj — kv cos 8) 

v 



(3.8) 
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In Ref. g it is stressed that the macroscopic description is meaningful only for distant 
collisions, and a framework which properly includes close collisions is formulated. Within 
such an approach the energy loss is UV convergent, but receives a contribution from the 
'hard' domain k S> T. In QCD the running of a s improves the UV convergence, and 
the stationary rate of energy loss (—dE/dx)^ = —AE oc (L)/L is expected to be actually 
dominated by the 'soft' k ~ m_D region We will come back to this point in the end of 
section 3.2. In the absence of a rigorous treatment with running coupling, we choose to 
focus on the difference between —AE(L) and the standard stationary law (|3.8| ), 

d(L) = —AE(L) + AE^L) (3.9) 
d(L) [ d 3 k ( I - cos {kL cos 6) 2 f°° dw 



C R a s J 2vr 2 1 k 2 + 



m 



+ v 2 [ \k 2 cos 2 9 p L + uj 2 s'm 2 9 p T ] 



sm 2 ((uj-kv cos6)L/(2v)) ttL 

2 - : ^= 6(u - kv cos 6 j 

[uj — kv cos 0) z v 



ind 

(3.10) 



which turns out to be UV convergent (see Appendix A). As discussed in the next section, we 
also expect the main feature of the function d(L) (which we evaluate with fixed coupling), 
namely its limiting value when L — > oo, to be unaffected by the running of a s . 

3.2 Numerical results and discussion 

In this section we first discuss the main features of d(L), in particular its large L limit, and 
the domain of validity of our calculation. (The mathematical properties of d(L) are studied 
in the Appendices.) Consequences on the phenomenology of the collisional loss —AE{L) 
are then presented. 

The large L limit of d{L), 

doo= lim d{L)= lim [—AE(L) + AE^L)] , (3.11) 



is calculated exactly in Appendix B (see Eqs. (B.13), (B.15) and ( |B.17| )) and can be accu- 
rately approximated by 

doc~-C K a s m D (l + v / 2(7-l) N ) . (3.12) 



Our central result is that doo scales as 7 = l/Vl — v 2 when v = p/E — > 1. The largeness 
(and negative sign) of d^ for large quark energies will be the main reason for the important 
delay of the stationary regime. The L-dependence of d(L) is presented in Fig. [I] for a fast 
charm quark. The function d(L) is found to be negative for all L, and Fig. |l]b shows the 
increase (in magnitude) of doo with the quark energy, as found analytically. 

The observed values of d^ can be translated to an estimate of the time scale charac- 
teristic of the transitory regime by defining the retardation time 

t ret = doo/ (dE/dx)oo ■ (3.13) 

For a quark with p = 10 GeV, -doo ~ 1 GeV and -(dE/dx)^ ~ 0.1 - 0.2 GeV/fm 
yields t rct ~ 5 — 10 fm. This number should be considered as a rough estimate, due to the 
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Figure 1: The function d(L) defined by (3J3) and ( 3.1 0| ) for a charm quark (M = 1.5 GeV) 
produced in a QGP of temperature T = 0.25 GeV, as a function of the distance L travelled in the 
plasma. We use a s = 0.2 and n/ = 2 in the expression of the Debye mass m 2 D = 4ira s T 2 (l +n//6), 
giving mo — 0.46 GeV. (a) d(L) for a quark momentum p — 10 GeV (full line). The contribution 
from initial bremsstrahlung (see section 4 and Appendix B) to d(L) is represented by the dashed 
line, (b) Dependence of d{L) on the charm quark momentum p. The thin straight lines give the 
values of d x for the corresponding values of p. 



numerous approximations used in our theoretical model. But we stress that the retardation 
time may be of the order of several fm for p = 10 GeV, and that it scales with the quark 
energy for large energies. 

The scaling in 7 of d^ (when 7 — > 00) results in a similar scaling of the retardation 
time, t re t ~ 7/ m D, and this has the following simple interpretation. The stationary regime 
for energy loss sets in when the regions of the plasma at a transverse distance ~ 1/mc, 
polarized by the parton current, start to retroact on the parton. When the parton velocity 
v is small, this takes a time t re t ~ l/m£>. When v — ► 1, the latter scenario can occur only 
if the parton is not ahead of the relevant polarized regions after the time i re t- At initial 
time t = 0, this requires the parton to polarize a domain which is ahead of it, at an angle 
6 ~ \/l — v 2 = I/7 1 with respect to the direction of v. The time t ret corresponds to 
the time necessary to send a signal at a transverse distance ~ 1/mu in the direction 8, 
leading to t ret ~ l/(mD0) for 9 <C 1, hence t TCt ~ ^/rri£, when v — > 1. 

We now discuss the domain of validity of our calculation. As shown in Appendix B 
(see ( B.18 )), the typical values of k contributing to cfco are k ~ O (jniD)- Since using the 
HTL gluon spectral densities in ( |3.10| ) a priori requires k <C T [l5|], our calculation 
of doo would seem to be justified, in the perturbative framework mo esc gT ^ T, only 
provided 7 is not too large, 7 <C 1/(7. However, as can be easily inferred from Appendix 
B, the dominant contribution to d^ arises from \uj 2 — k 2 \ ~ m? D (with uj ~ k ~ r yni£)), i.e. 
from the region where K = (w, k) is close to its mass-shell. It is known ]i~9| , |20| ] that in this 
region of low virtualities the HTL gluon propagator is a very good approximation to the 
exact propagator even in the domain k ^> T. Thus our calculation of d^ might be justified 
for all values of 7. 

Contrary to the position of the asymptote d^,, the small L behaviour of d(L) seen on 
Fig. H should not be physically sound. The small L behaviour of d(L) can be obtained 
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Figure 2: Collisional energy loss of a charm quark produced in a QGP as a function of the distance 
L travelled in the plasma, for p = 10 GeV and p = 20 GeV. The total result (full line) is compared 
to the linear law (— AEoc) corresponding to ( 3.14 ), and which would arise in Ref. [Q by evaluating 
the slope to logarithmic accuracy. As in Fig. [l]a the dashed line singles out the contribution to 
—AE(L) from initial bremsstrahlung (see section 4). 



for instance from ( p.4[ ), and is of the form d(L <C m^ 1 ) oc a 2 T 2 L\og (itidL), where the 
logarithm arises from an integral ~ dk/k. Thus the small L limit is sensitive to large 
virtualities \uj 2 — k 2 \ ~ k 2 ^> T 2 , where the HTL approximation becomes invalid. 

In order to qualitatively illustrate how d(L) delays the collisional energy loss, we have 
to add the stationary contribution — A£' 00 (L) given in (3^8). In the end of section 3.1 we 



pointed out that (|3.8|) is ill-defined as it stands, because of a logarithmic UV divergence. 
We mentioned that this divergence would be absent, either with a proper treatment of close 
collisions with fixed a s [f|, or due to the running of a s ||, both approaches differing by the 
hardness of k contributing to — (dE/dx)^. This stresses the need for a rigorous treatment 
with running a s . In the absence of such a treatment, and since we need the stationary 
loss — A£' 00 (L) only for a qualitative illustration of the retardation effect, we will use the 
'standard' result derived in Ref. Q, which reads to logarithmic accuracy (which is sufficient 
for our purpose): 



dE\ 

dx) { 



-AE^L) _ C R a s m 2 D 



L 



where k n 



1 

v 

Min 



1 



2v< 



log 



1 + v 



log 



A* T , 




(3.14) 



(3.15) 



and m g = m£)/V3 is the gluon thermal mass. 



Our results for —AE(L) given by (|J), (|3T0| ), (|1J) and (|1|) are shown in Fig. |, 
illustrating the main qualitative feature, i.e. the important delay of the stationary regime. 
The retardation time (3.13), which corresponds to the value of L where the asymptote of 
—AE(L) cuts the horizontal axis, is close to the intersection of the curve — AE{L) with 
this axis. 

As another consequence of the large negative values of d(L), we observe on Fig. | that 
the induced energy loss — AE(L) is negative at relatively small L, and thus corresponds to 
an (induced) energy gain. 
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Two effects may explain the latter observation, which goes beyond our initial expecta- 
tion of a delayed but however positive energy loss. First, we recall that current conservation 
requires the energetic parton to be produced with at least another parton. In the simple 
case of a color singlet dipole we considered (see (^q)), we expect the dipole to separate 
more easily in the medium than in vacuum, due to charge screening. Secondly, the current 
created at t = produces radiation 4 . As discussed in detail in section 4, in the medium 
the radiated energy corresponds to the excitation of the QGP plasmon modes and is re- 
duced compared to vacuum (the induced radiated energy is negative) , as might have been 
expected for massive modes. However, the retardation effect is only partly due to this dif- 
ference between in-medium and vacuum radiation, as seen in Fig. |l|a and further discussed 
in section 4. 

We end this section by discussing which features of our calculation might be affected 
by the running of a s . 

First, the standard result 0] for — A£ ; 00 (L) we have used arises from the logarithmic 
interval m g <C k <C fe max (see ( |3.14| ) and ( |3.15| )). In QCD, with a running coupling 
evaluated at a scale on the order of the gluon virtuality, we expect {—dE/dx)^ to be 
oc J dk a s (k) 2 /k oc J dkj(k log 2 A;), showing that (—dE/dx)^ depends negligibly on E 
when E — >• oo || . (In this sense the log k max oc log E factor in ( p. 14 ) is an artefact of the 



fixed coupling approximation.) Most importantly, (—dE/dx)^ is actually dominated by 
the soft (infrared) region k ~ mo when E 3> ran ||, and we expect the calculation of 
( — dE/dx)^ within a macroscopic description and with running a s to be self-consistent. 

Secondly, the running of a s should affect the behaviour of d(L) at small L < l/mo- We 
discussed previously that since the slope of d(L) at small L is of the form ~ a 2 dk/k, 
it cannot be consistently derived within a macroscopic description. This should not be the 
case with running coupling, since we expect the slope to be rather ~ X^f dk/(k log 2 k), 
which is dominated by k ~ mo and independent of L when L — > 0. Thus the sharp be- 
haviour d(L — > 0) oc LlogL obtained for fixed a s should become d(L — > 0) oc —a s (m£) )T 2 L 
in the case of running a s . As a consequence, the energy loss —AE(L) = — AS 00 (L) + d(L) 
shown in Fig. [2| should be at most 5 ~ O (L) at small L. Thus, the relatively important 
induced energy gain at small L seen in Fig. [2] might be strongly affected by the running of 
the coupling. In contradistinction, we stress that the large L limit d^ of d(L) should not 
be affected by the running of a s , since doo depends on low virtualities \lu 2 — k 2 \ ~ m 2 D . 

4. Induced radiation 

As already mentioned, and as is well-known, the sudden acceleration of the energetic parton 
at time t = comes along with bremsstrahlung radiation. The quantity —AE(L) which 
we called induced collisional energy loss actually contains this contribution. In this section 



4 The induced radiated energy is part of —AE(L) as we defined it (see section 4). 

5 In fact, for running coupling the slope of —AE(L) is expected to vanish when L — » 0. This can be seen 



by expanding the bracket in (2.9) to order L. The UV convergence of the integral over k being ensured 
by the running of the coupling, the integral over w can be performed by closing the contour in the upper 
half-plane where there is no singularity and identically vanishes. 
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6 (rad) 



Figure 3: In-medium angular radiation spectrum (broad full line) ( |4.2| ) for L = 5fm, p = lOGeV, 
k — 0.5 GeV (left) and k = 1 GeV (right). The thin full line represents the in-medium spectrum in 
the L — > oo limit (4.2). The broad and thin dashed lines represent the same spectra in the vacuum, 



given respectively by (4^4) and (|4.q). 



we single out this radiative component from —AE(L) and emphasize that it plays only a 
finite part in the retardation effect studied in section 3. 

In order to obtain the contribution W(L) to —AE(L) originating from radiation, we 
single out in (|3.4|) the contribution of the plasmon poles. This is done by using the expres- 



sion of the gluon propagators close to their poles [15], i.e. at or ~ i^i(k), 



A s (u, k) 



-z s (k) 



(u + irj) 2 — 0J 2 (k) 



-z s (k) 



a* (A:) 



in sgn(a;) 5(uj 2 — uj 2 (k)) 



(4.1) 



and picking only the <5-term in the latter. From ( |3.4f ) we obtain the total (i.e., we do not 
subtract the vacuum contribution from it) in-medium radiation spectrum: 

k 2 ^2, 



dW(L) C R a s 



dk d cos l 



Crc< s 



IT 



z L (k) 



2 sm 2 ((oj L {k) - kv cos 9) L/(2v)) 
COS a - 



2tt 



z L (k) 



k 2 



+ z T (k) sin 2 9 



(cos 9 -uj L {k)/{kv)f 
sm 2 {{uj T (k) - kv cos 9) L/(2v)) \ 



(cos 9 - LU T (k)/(kv)Y 



7 



cos 



(cos I 



z T (k) 



sm 



(4.2) 
(4.3) 

5 fm, 



J L \rvj poi/ — u) L(k) / (kv)) 2 ' ± v ' (cos 9 — ujx(k) I (kv)) 2 

where the functions u)L y x(k) and 2l,t(&) can be found in Ref. [jig] . 

The differential spectrum in the polar angle is represented in Fig. || for L 
p = 10 GeV, and for two values oik, k = 0.5 GeV and k = 1 GeV. The vacuum contribution 
is obtained by setting zl = 0, zt = 1, ^h{k) = wy(fc) = in (^), yielding: 

C^cvs . 2 n sin 2 ((A; - kv cos 9) L / (2v)) 



dW(L) 



dk d cos 9 



sm 



7T 



(cos 9 -l/vf 



C R a s 



v 2 sin 2 i 



L^oo 2-7T (1 — VCOs9) x 



(4.4) 
(4.5) 



As expected, (4J3) corresponds to the vacuum bremsstrahlung spectrum of a single charge 
suddenly accelerated at t = 0. We remark that the latter spectrum contributes to the 
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radiative energy loss induced by parton multiple scattering |6|, 0, §, as the zeroth order 
in the number of parton rescatterings ('self-quenching' j8|). As already mentioned, the 
corresponding in-medium radiation spectrum ( f4.3[ ) arises from medium polarization, and 
thus differs from the spectrum induced by (at least one) rescattering calculated in |], ^, [8| . 
The (large L) in-medium radiation spectrum ( |4.3| ) was previously obtained in Ref. |2l| 
within a diagrammatic approach. 



In vacuum the spectrum ([Oj) is modified by parton cascade and hadronization, but 
the angular pattern fl4.5| ) in principle affects the distribution of final state hadrons. In a 
QGP of very large size, this distribution will rather be sensitive to the large L in-medium 
spectrum (|4.3|). Similarly, the modification of the finite L diffraction pattern when one goes 



from vacuum (spectrum (4.4)) to the medium (spectrum ( |4.2| )) may be of some relevance 
to discuss particle production induced by a fast parton travelling a distance L in the QGP. 
When k is large, the spectrum ( [4.2|) approaches the vacuum spectrum ( |4.4| ). This can 
already be seen for k = 1 GeV (Fig. || right), although some attenuation is still visible at 
small 9. When k decreases (Fig. || left), a strong distorsion shows up as ( |4,2| ) becomes 
sensitive to the plasmon modes of the QGP. 



However the spectrum (|4.2j ) is not directly observable (not speaking of parton cascade 
nor hadronization). In a realistic situation, the parton propagates in the thermal medium 
for < t < L/v, and in vacuum for t > L/v. Thus the full picture should incorporate 
both stages. In addition, the fast parton produces transition radiation when going through 
the discontinuity between the thermal medium (ei,T 7^ 1) and the vacuum (cl,t = !)• 
This transition radiation must also be taken into account in order to obtain the complete 
radiation spectrum. We postpone to a future work the study of transition radiation in our 
context, as well as its influence on the overall energy loss. In particular, it is too early 
to possibly relate the angular spectrum ( |4.2| ) (see Fig. ||) to the humpbacked azimuthal 
distribution of particles produced back to the leading jet observed at RHIC p2| . However, 
the angular radiation spectrum (^^) shown in Fig. ^ suggests that (at least some component 
of) the particle yield in heavy ion collisions should exhibit a diffraction pattern arising from 
the finite length of the medium along which collisional energy loss occurs. 

Finally, we stress that in our calculation the radiated energy (obtained by integrating 
( |4.2| ) over the radiated gluon energy and emission angle) is smaller in the medium than 
in vacuum. This feature could be foreseen by comparing the in-medium and vacuum 
differential spectra in Fig. [3|. Thus the contribution of induced radiation to what we 
defined as the induced collisional energy loss — AE(L) is negative. Numerically, we find 
that the bremsstrahlung contribution always accounts for less than 50% of the difference 
d(L) between —AE(L) and the asymptotic result — AE^ (see Figs. |l]a and |J). In the 
end of Appendix B we show that the relative contribution from bremsstrahlung to d(L) 
reaches exactly 1/2 when L — > oo and in the ultrarelativistic limit v — > 1. This shows 
that the retardation effect studied in section 3 cannot be solely attributed to the radiative 
component of —AE{L). 
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5. Conclusion 



We have studied the collisional energy loss —AE(L), in the fixed coupling approximation, 
of an energetic parton travelling the distance L in a QGP, and initially produced (at 
t = 0) in the medium. Compared to previous estimates which assumed the parton to be 
produced at t = — oo, —AE(L) is strongly suppressed up to L ~ t ret , where the retardation 
time t ret scales with the parton momentum p at large p. For p = 10 GeV we roughly 
estimated t ret ~ 5 — 10 fm. We stress that this estimate is only indicative due to the 
various approximations used in our theoretical model, in particular the small coupling 
limit g < 1. Also, the running of a s is expected to affect both the small L behaviour 
of —AE(L), and the asymptotic stationary rate (—dE/dx)^, as discussed in the end of 
section 3.2. A rigorous treatment with running a s would be needed to obtain a better 
quantitative estimate of i re t • 

We believe our results could be relevant to jet quenching phenomenology, since it 
has recently been argued || [ll| that collisional energy loss should be reconsidered as an 
important source of energy loss. In addition to the suppression of —AE(L), we find that the 
asymptotic linear behaviour of —AE(L) is delayed to about L ~ i re t. In particular, using 
a stationary energy loss rate —dE/dx makes sense only for quite large values of L > t Te t- 

The suppression and retardation of —AE(L) are encoded in the difference d(L) = 
—AE(L) + AE 00 {L) between —AE(L) and the 'standard' stationary result —AE^L) oc L. 
As we have shown, d(L) is a well-defined (UV convergent) quantity, and our main result - 
the largeness of the retardation time - is independent of the uncertainties on (—dE/dx)^ 
and on the precise shape of —AE(L) at small L. The main reason for the large magnitude 
of the retardation time is the scaling d^ cx —^m£> when 7 = E/M — > 00, where doo is the 
large L limit of d(L). We showed that the (induced) bremsstrahlung arising from the hard 
parton being suddenly accelerated at t = 0, formally included in d(L), contributes only 
partly to d^. Thus the retardation effect cannot be solely attributed to initial radiation. 
The scaling in 7 of doo (when 7 — > 00) results in a similar scaling of the retardation time, 
^ret ~ l/ m Di an d a physical interpretation of this fact is given in section 3.2. 



Finally, we mention that in order to explain the observed dependence p3| in azimuthal 
angle <ft (with respect to the reaction plane in a heavy ion AA collision) of the nuclear 
modification factor Raa, taking into account the geometry of the collision is not sufficient 
1 24]. A length scale L ~ 2fm has to be introduced, below which the parent parton of the 



high pt jet or hadron is assumed to be insensitive to energy loss. In Ref. [24] this parameter 
is interpreted as the formation time of the plasma, but it is also stressed there that 2 fm is 
quite large compared to the values ~ 0.2 fm usually taken for the plasma formation time. 
Our study suggests that this parameter might instead hint to the possibility of a negative 
loss — AE co u(L) before the stationary regime, partially compensating the radiative energy 
loss — AE rac [(L) induced by rescattering of the energetic parton. 
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APPENDIX 



A. Ultraviolet convergence of the function d(L) 

Here we show that d(L) given by ( 3,10| ) is ultraviolet convergent. For this purpose, it is 
sufficient to prove that the angular integral 



c?cos 8 / (cos 9) 



k 2 (x — vcos8) 2 kv 2 \v 



(A.l) 



where x = oj/k and /(cos 9) = cos 2 9 or sin 2 9, is of order ~ O {l/k 2 \ when k — » oo. Indeed, 
the second term of ( |A.1[ ) is of order ~ O (1/&), leading to the logarithmic UV divergence 
of — AE' 00 (L) (see d3~8[)). With the change of variable 



u = a (x — v cos 9) ; a = — 

2v 



(A.2) 



we obtain 



L 



kv 2 

where we define 



a(x+v) 



du- 



sm u „ (x u 



a(x—v) 



v av 



ft ^^(v-lxDf - 



X 



(A.3) 



h 



L 

kv 2 



a(x+v) sin 2 u 

du — k — 

a(x— v) 



V- 



tt9(v — \x\) 



ra(x+v) 



kv 2 



du- 



sm 2 u 



a(x—v) 



11- 



f 



X U \ 

v av) 



(A.4) 
(A.5) 



The behaviour of I\ when a = kL/(2v) — > oo is easily found by treating separately 
the two cases |x| > v and |x| < v. In the first case we can replace sin 2 u — > 1/2 in the 
integrand and we find 



q^oo kv 2 \v/ax 2 — v 2 



1 



V J x 2 — V 2 



(A.6) 



which is ~ O (l//c 2 ). In the second case, |x| < v, we write 



*(x+v) gin 2 u 



du- 



ii 



a(x—v) 



du- 



sm 2 u 



a(x+v) 



*{x-v) s[n 2 u 

du — k — 



(A.7) 
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where the replacement sin 2 u — > 1/2 can be made in the r.h.s. when a — > oo, leading again 
to ( |A,6| ). Note that when a — > oo, the expression (A. 4) is thus equivalent to 



L „ (x\ /■«(*+«) , „ / 1 

' a(x— v) 



In the ^2 integral defined in (|A.5| ), the contribution from u <C a is negligible, showing 
that typically u ~ a. Thus we can replace sin 2 u — > 1/2, yielding: 

„ ~ 'f^pWf/p.i)./^ . (A . 9) 

a-*oo 2kv z J a r x -v) \u z / \ \v avJ \v 



We infer from the above that the behaviour of (A..1) when a = kL/(2v) — > oo is 
obtained by the following replacement: 



k 2 (x — v cos 9) 2 kv 2 \v 



a— >oo 



p l, gF^ ' ' <A ' 10) 



The r.h.s. is ~ O (l/fc 2 ) , to be compared to O (L/k) (second term of ( |A.1| )). This completes 
our proof of the UV convergence of the function d(L). 

B. Large L limit of d(L) 

Here we evaluate the limiting value of d(L) (defined by (|3.9|)) when L — > oo, 



doo= lim d(L) = lim [-AE(L) + AEoo(L)] . (B.l) 
B.l Total result for doo 

Using (|3.4| ) (with the bracket expressed as in fl2.10| )) we obtain: 



d(L) „ f d 3 k 



x (,«v.- ! ™-.^Mp(i]|. M 



We now use the fact that if we replace L by — L in (^4), the integral over u, performed 
by closing the integration contour in the upper half-plane, vanishes identically because the 
singularities in u all lie below the real axis: 

d 3 k 



o = iv 2 f il f°° \k 2 C os 2 e& L {u,k)+uj 2 sm 2 e& T {uM 

J 4vr d y_ oc to 



Lad 



-, Si K, Wv + 2 - i g^Mp | i . ,B,, 



Adding (JBJ) to ([B.2D we get: 



CrOs 



-2i** KyL '%(J-)\ . (B.4) 
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-S'(x), we obtain in the 



Using the first of the identities ( 2.11] ), and then 5(x)P(l/x) 
L —* oo limit: 

VRe [t^f C ° d ^- [k 2 cos 2 A L (co, k) + u 2 sin 2 9 A T (co,k)]. nd 5'(K.V) . 

(B.5) 



doo 
Cro s 



2tt 3 



In <J'(if.y) we trade the derivative d/duj for d/dv: 

doo 2 ^ 1 „ /" d 3 fc 

7T?; 



Using (|3.3|) we have 



- — Re / — 5- \A L (kvcos9,k) + v 2 sin 2 9 A T (kv cos9,k)]. 
ovv J 2tt 6 1 J 1 

n L (x) 



[Ai(fei,fc)]i 
[AT(^x,A;)] ind 



" fc 2 [A; 2 + n L (x)] 

-n T (x) 



£; 2 (x 2 -l)[P(x 2 -l)-n T (x)] 
and denoting x = v cos 6 the equation ( |B.6[) becomes: 



(B.6) 

(B.7) 
(B.8) 



<K odl f v dx rdk Re i n z (aO t> 2 -x 2 n T (x) 



<9v 



U jo 



7T 



A; 2 + U L (x) 1 - x 2 k 2 (l - x 2 ) + n T (x) 



We finally perform the integral over k, 



d d 1 f v 

rriDV 2 — ^ I dxRe i \/Ul(x) 



Cro s 



dv v 2 Jq 



where 



U L (x) = U L (x)/m 2 D = 1 - - log 



v 2 — x 2 



(1 - x 2 fl 2 



x + 1 



U T (x) 



(B.9) 
(B.10) 



x- 1 



n T (x) = U T (x)/m 2 D 



x 

Y 



x 2 -l 
2x 



log 



x + 1 
x-1 



The longitudinal contribution to (B.10) reads: 

dooL = -C R a s m D A L (^) 



(B.ll) 



(B.12) 



A ^) = \J dxReWl-|log 



x + 1 



x — 1 



— -ReJl--log 



v + 1 



v-1 



We check numerically that A/, (7) is a smooth increasing function of the parton Lorentz 
factor 7 = — v 2 , increasing very slowly above 7 = 10, and saturating when 7 — » 00: 

= 1, Ai(10) ~ 1.3, A L (1000) ~ 1.5, A L (oo) ~ 1.814. 
For the transverse contribution to ( B.10| ) we find: 

dooT = -C R a s m D At (7) ( B - 14 ) 



2 

A r ( 7 ) = - / dx— 



!)3/2 



ReW — 



x 2 -l 



2x 



log 



x+1 



x-1 



+ 



mx{x 2 — 1) 



(B.15) 
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Figure 4: The exact function — (1^/(0 Ra s m£,) (full line) and its approximation 1 + \/2(7 — 1) 



(dashed line), as a function of 7 = E/M, see (B.17) 



When v — * 1, a singularity at x — ► 1 appears in the integrand of ( |B.15| ), and we easily 
derive the asymptotic behaviour ^4t(t) — V^T when 7 — > 00, giving: 

dooT ^ -\/2C R a s m D ^ . (B.16) 

U— »1 

Numerically, -At (7) (which satisfies = 0) is quite well approximated by the linear 

form v2(7 — 1) for all values of 7 > 2, with an accuracy of 25% for 7 = 2 and improving 
for increasing 7. 

Adding flRl2| ) and pUD we obtain 



doo = -C R a s m D {A L {>y) + A T ( 7 )) ~ -C R a„m D 1 + V2( 7 - 1) , (B.17) 



where the latter approximation has the correct limits at 7 — ► 1 and 7 — > 00, and can be 
checked numerically to be accurate to better than 12% for all values of 7 (see Fig. ||). 

It is instructive to find what are the typical values of k and 9 which contribute to 
the leading behaviour ( |B.16| ) of d^. In (|B.9| ) the transverse contribution is dominated by 
x = v cos 9 — > 1 when u — ► 1, i.e. 9 < 1. Using Hr(x — ► 1) = m 2 D /2, the leading part of 



(B.9) can be written as 



a -r f 1 d ° 2 f°° dk m D mm 

Ras Jo ez + l/jzJo n k 2 (9 2 + l/ 7 2 )+m 2 D /2 ■ ^-i«J 

The latter is easily checked to arise from 9 ~ I/7 and ~ 7m_D (as well as to yield the 

result ( pnep ). 

B.2 Contribution from radiation to doo, m the v — > 1 limit 

In order to see how the bremsstrahlung induced by the initial acceleration of the parton 
affects the retardation effect, we single out in d^, given by ( |B.17|) , the contribution from 
radiation. This contribution is obtained by subtracting ( |4.5| ) from (|4.3| ), and then 
integrating over k and 9. Here we concentrate on the transverse contribution Wqot which 
turns out to be dominant in the v — > 1 limit, 

Woot f°°dk f 1 f z T (k) sin 2 9 sin 2 9 \ 

~nZ / dcos9\ „,,w,,„u 2 - 7—7 ; /„v 2 > • (B.19) 



Cija s 7 2tt J_x \(cos9 - u; T (k)/(kv)) 2 (cosO-l/vf 
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When v — > 1, the above integral is dominated by the domain 6 <C 1, k S> mp. For k ^> 
we have uoT{k) ~ k+m^/ (2k), where moo = mo/vl is the asymptotic gluon thermal mass, 



and zx(k) ~ 1 |L5|] can also be consistently used in (|B,19 ). Approximating the integrand 
in ( B.19D we obtain 



dd 2 e 2 =-^: — ^ , • ( B - 2 °) 



C B a s t^l ir J k* J Q ( 02+ i \ 2 ^2 + i + 



In the ultrarelativistic 7 — ► 00 limit, the typical values of k and 6 in the latter integral 
are k ~ 777(00 and ~ I/7. The calculation is now straightforward and we obtain for the 
leading term: 

Woot ^ -C R a s m Dl /V2 . (B.21) 

Xhis is exactly half of the full result for d^, see ( gig ) and (|B. 17| ). 

For completeness we quote the result for the longitudinal contribution to Woo, 

tit 2CRa s m D 3/2 

W^ooL ~ (log 7) 7 , (B.22) 

U-+1 07T 



which is indeed subleading compared to (B.21) when v 
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